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Abstract 

The proper definition of entropy is fundamental to the relationship between statistical mechanics 
and thermodynamics. It also plays a major role in the recent debate about the validity of the 
concept of negative temperature. In this paper, I analyze and calculate the thermodynamic entropy 
for large, but finite quantum mechanical systems. A special feature of this analysis is that the 
thermodynamic energy of a quantum system is shown to be a continuous variable, rather than 
being associated with discrete energy eigenvalues. Calculations of the entropy as a function of 
energy can be carried out with a Legendre transform of thermodynamic potentials obtained from a 
canonical ensemble. The resultant expressions for the entropy are also able to describe equilibrium 
between quantum systems having incommensurate energy-level spacings. This definition of entropy 
preserves all required thermodynamic properties, including satisfaction of all postulates and laws 
of thermodynamics. It also demonstrates the consistency of the concept of negative temperature 
with the principles of thermodynamics. 
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I. INTRODUCTION 


In recent years, there has been a renewal of interest in the fundamental dehnition of 
the thermodynamic entropy in statistical mechanics. An important source of that interest 
has been the controversial claims that the Gibbs “volume” entropy (dehned below) is the 
only valid de£nition[THn], and that the concept of negative temperature is inconsistent with 
thermodynamics [711T5] ■ Earlier work on classical systems has shown how a valid dehnition of 
entropy is consistent with negative temperatures, giving correct thermodynamic predictions 
when the volume entropy does not [T61 - I23] . 

The entropy of quantum systems is complicated by a discrete energy spectrum, which 
has led many to regard the entropy as being dehned only at energies corresponding to 
eigenstates [211 121], or as constants between eigenvalue energies (step functions) jH-fO] . In 
this paper, I argue that the thermodynamic energy of macroscopic quantum systems has a 
continuous spectrum of values, even though the eigenvalues of the energy form a discrete set. 
To briehy summarize my arguments, macroscopic measurements never put a many-body 
system into an energy eigenstate. The small interactions that are necessary to establish 
thermodynamic equilibrium between two or more systems are sufficient to alter the discrete 
energy spectra of the isolated systems. This means that the projection of an eigenstate 
of the composite system of interacting subsystems onto the energy levels of an isolated 
subsystem will never produce an eigenstate of the isolated subsystem. An isolated system 
will necessarily be described by an ensemble of quantum states, including contributions from 
many energy eigenvalues. 

I begin by demonstrating in Section |TT| that the thermodynamic energy of an isolated 
quantum system that was previously in thermal contact with another system is a continuous 
variable. Section [ITT| reviews the canonical probability distribution of energy in a system that 
has been in contact with a thermal reservoir, that is, with a very large (formally, inhnite) 
system. I analyze the canonical distribution in terms of a Massieu function |2Bl 127] . in order 
to be able to carry out calculations for either a monotonic or a non-monotonic energy density 
of states. The essential features of an ensemble for an isolated system are then shown to 
depend on an internal temperature variable, which leads to a method for calculating the 
entropy of a macroscopic quantum system. As examples, explicit expressions are derived for 
the entropy of quantum simple harmonic oscillators and two-level objects. The entropy of a 
system of two-level objects has a non-monotonic density of states, and serves to demonstrate 
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that negative temperatures are consistent with the principles of both quantum statistical 
mechanics and thermodynamics, as previously shown for the classical case in Ref. |231 


II. THERMODYNAMIC ENERGY IN QUANTUM SYSTEMS 

To see the unique way in which interactions affect thermal ensembles in quantum statisti¬ 
cal mechanics and lead to the thermodynamic energy being a continuous function, consider a 
composite system with M macroscopic (many-body) subsystems. The macroscopic observ¬ 
ables for the composite system are just the set of 3M variables {E,V,N} = {Ej, Vj, Nj\j = 
1, ■ ■ ■ , M}, where E,V, and N denote the sets of energy, volume, and particle numbers. To 
reduce the proliferation of subscripts, the equations are written for a single type of parti¬ 
cle, but the generalization to many types of particles is obvious. For quantum spins, the 
magnetization will also be included in the macroscopic variables. 

Each macroscopic subsystem j will be characterized by a Hamiltonian Hj, when it is iso¬ 
lated from other subsystems. To establish thermal equilibrium between subsystems, direct 
interactions between particles in different subsystems are necessary. These additional inter¬ 
actions will be denoted by the interaction Hamiltonians L = {Lj^k\j = 1, ■■■, M — 1] k = 
j,-- - ,M;}. These interactions between macroscopic subsystems prove to be essential in 
deriving the quantum entropy. 

The extension of the dehnition of entropy to quantum systems follows the basic pattern 
of considering thermodynamic interactions between subsystems of a composite system that 
I discussed for classical composite systems in earlier work|21j. However, there are some 
subtleties associated with the discreteness of the spectrum of energy eigenvalues for hnite 
systems. In particular, the effects of interactions between subsystems, denoted by L, must 
be treated carefully. I begin by introducing a fairly standard notation for quantum states 
and equilibrium quantum ensembles, hrst for isolated subsystems, and then for interacting 
subsystems. 


II. 1. Isolated quantum subsystems 

First consider the case in which L = 0, so that every macroscopic subsystem is isolated. 
Denote the quantum Hamiltonian for the j-th subsystem by Hj, and the corresponding 
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quantum numbers by Vj. The energy eigenvalue equation is then 


Hj\pj) = ( 1 ) 

where \i'j) denotes an eigenstate. I assume that for each subsystem j, the set of \i'j) form a 
complete, ortho normal set. 

A general quantum state of subsystem j will be denoted by It can be expanded in 
the eigenstates of the subsystem, 

\'^j) = (2) 

where the Cjy^.’s are complex numbers, Cy. = \cuj \ exp 

Cuj = ( 3 ) 

If l^pj) is properly normalized, we have 

Assuming that there is a probability distribution, Pj{'il!j), dehned on the set of quantum 
states, then the average of an arbitrary operator A will be given by 

hi 

where the double angular brackets indicate that both a quantum expectation value and 
an ensemble average are being calculated, and the summation symbol is intended to cover 
whatever sums and integrals are appropriate. Using Eq. (|^, ((A)) can be written as 

«'4)> (6) 

{Cj} H 

where {cj} indicates a set of values of c^. corresponding to the quantum state \'ipj), and the 
sum over {cj} is the sum over all quantum states in the ensemble. 

In equilibrium, all phases, are equally probable. Integrating over the phases gives a 
factor of zero for all cross terms, leaving the expression 

«-4)) = E E (7) 

’^3 {Cjl 
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This expression can be simplified by writing 

~ I ) ( 8 ) 

{cj} 

so that Pivj) > 0. This gives 

{{-4)) = 5^P,(>'i){>'yl^l>',). (9) 

^3 

Further noting that ((!)) = 1, we see that 

^3 

SO that Pj{i'j) acts very much like a probability. Averages over the equilibrium quantum 
ensemble can be computed as if Pj{vj) were a probability. It is quite common to speak 
loosely of Pjii'j) as being the probability of the subsystem being in the |z/j) eigenstate, but 
it is well known this is not correct. Pj{vj) is the probability that a precise measurement of 
the energy would put the subsystem into the |z/j) eigenstate. Since a macroscopic subsystem 
is generally in a state of the form given in Eq. ([^, and macroscopic measurements cannot 
put a subsystem in an eigenstate, the true probability of a real macroscopic subsystem being 
in an energy eigenstate is zero. This leads to the question of what we do know about the 
distribution of quantum states in interacting subsystems. 

II.2. Interacting qnantum snbsystems 

An important feature of interacting quantum subsystems is that the spectrum of energy 
eigenvalues is not related to the corresponding spectra for the isolated subsystems in a simple 
way. The composite subsystem has new eigenfunctions with new eigenvalues that depend 
on the interactions L 7 ^ 0. As an example, consider two interacting subsystems with the 
following Hamiltonian. 

Pi ,2 = Hi + H 2 + Li2 ( 11 ) 

The eigenvalue equation for this Hamiltonian can be written as 

Hl,2\l^l,2) = ei,2(i^l,2)|^^l,2), (12) 

where |i^i, 2 ) denotes the eigenstate, and ^ 1^2 the quantum number. The eigenstates \i'i, 2 ) 
form a complete, orthonormal set. 
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A general quantum state of the composite system of interacting subsystems will be de¬ 
noted by Id'ig). It can be expanded in the eigenstates of the composite system. 


l '\,2 

where the C'l.^j’s are complex numbers. 




If Id' 1 , 2 ) is properly normalized, we have 

^'1.2 


(13) 


(14) 


(15) 


II.3. Energy distribution of an individual subsystem due to interactions 

Given the eigenstates and eigenvalues of the composite system of interacting subsystems 
and the probability distribution ^ 1 , 2 (^ 212 ), we can compute the energy distribution in the 
individual subsystems. If the composite system is in a state ld' 12 ), the projection onto the 
eigenstate |z/i) of the isolated subsystem 1 is 

(Z2l|d'l,2) = 5^G^l,2(i2l|l2l,2) (16) 

1 ^ 1,2 

The expectation of the energy Ei of subsystem 1, given that the composite system is in 
state ld/ 1 ^ 2 ), is 




(17) 


Ml 


Because |i 2 i) is an eigenfunction of Hi, Eq. (01 becomes 


or, using Eq. (16), 


^ 1 ) (^1-^1 ^l)(^l ^ 1 , 2 ) 

(18) 

{Bi)»„ = ^|a„pi(>'iki,2)tBi(i^i) 

(19) 


I2l 


Averaging over the probability distribution of the states of the composite system gives 

((-^ 1 )) = X^-Pi,2(^2i, 2) i^l,2)P (20) 
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This can be simplified by defining 

Pl(^l) = -fl,2(^l,2) |C'j/i,2| \{^l\^l,2)?■ (21) 

1 ^ 1,2 

where 1 > Piivi) > 0, and Yhy^Viiyi) = 1- This gives 

({£.)) = ( 22 ) 

The value {{Ei)), including both quantum expectation values and ensemble averages, cor¬ 
responds to the thermodynamic energy, Ui, so that we can also write 




(23) 


The most important feature of Eqs. (22) and (23) is that they demonstrate that subsys¬ 
tems that interact within a composite system are not in energy eigenstates. This, in turn, 
implies that the thermodynamic energy of a quantum subsystem in equilibrium with another 
quantum subsystem can vary continuously, even though the energy spectra are discrete. Ui 
remains a continuous variable when subsystems in thermal equilibrium are separated. 

It is appropriate to compare the magnitudes of typical energy spacings with typical ther¬ 
mal fluctuations of the energy. In an example presented below, we will consider a subsystem 
composed of N two-level quantum objects, with a level spacing of e. Characteristic average 
energies are of the order of Ne, so that the relative size of the level spacing in comparison 
with the total energy is of the order of 1/N. In contrast, relative fluctuations of the energy 
are typically of the order of l/\fN. For a macroscopic subsystem with N ^ 10^°, the thermal 
fluctuations are roughly 10^° larger than the energy level spacing. 

It is also important that typical energies associated with the interactions L are expected 
to scale with the surface area of the subsystem, so that their relative magnitude goes as 
^-1/3 This is can easily be larger than the thermal fluctuations, and much larger than 
the energy-level spacing. As a result, there is nothing to prohibit quantum subsystems 
with incommensurate energy-level spacings from exchanging energy and being in thermal 
equilibrium with each other. Indeed, the explicit forms for the entropies of model subsystems 


that are derived in Section [IV] can be used to predict thermal behavior due to interactions 
between macroscopic subsystems with incommensurate energy-level spacings. 


Another consequence of Eq. (23) is that the concept of a microcanonical ensemble must 
be modified in quantum statistical mechanics, as discussed in the next subsection. 
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II.4. Separated quantum subsystems 


When two or more interacting classical subsystems are separated and isolated, they will 
each go into a microscopic state with a unique energy. Naturally, the value of that energy 
can only be predicted to within the thermal fluctuations of the energy before separation. 
Nevertheless, it clearly makes sense to describe an isolated classical subsystem by a micro- 
canonical ensemble with an exact value for the energy, regardless of its history. 

The situation when quantum subsystems are separated is quite different. Each subsystem 
will go into a quantum state that can be expressed as a linear combination of eigenstates, 
as shown in Eq. (|^. However, Eq. (23) shows this quantum state will not be in an energy 
eigenstate, so the subsystems cannot be described by the usual microcanonical ensemble 
with an exact value for the energy. As demonstrated in the next section, it is still possible 
to compute the fundamental relation, Sj = Sj{Uj,Vj, Nj), for each subsystem by using the 
canonical ensemble, which is valid for quantum statistical mechanics. Explicit examples of 


such calculations are then given in Section IV 


III. CALCULATING THE ENTROPY THROUGH THE CANONICAL ENSEM¬ 
BLE 


To determine the average thermodynamic energy we must calculate the “probabilities,’ 


that were dehned in Eq. (21). The simplest way to do this is through the canonical 
ensemble, which will give us a thermodynamic potential that contains all thermodynamic 
information for the subsystem of interest. 

A Legendre transform will then give the entropy as a function of the thermodynamic 
energy. Because I am interested in using a formalism that is also valid for subsystems with 
non-monotonic entropies, I will use a Massieu function to derive the entropy [2^ I7F] . Since 
a Massieu function might be somewhat less familiar than the more common Helmholtz free 
energy, it will be reviewed in the subsection following that on the canonical ensemble. 


HI.l. Canonical ensemble 

If we put the subsystem of interest, which we will pick to be j = 1, into equilibrium with 
a reservoir that has a known continuous density of states, ujr{Er), the value of piivi) will 





be determined by the proportionality: 


Pi(z/i) oc wh {Et - Ei{i'i )), 


(24) 


where Et = Er + Ei{vi) is the total energy of the reservoir and subsystem 1. Taking the 
logarithm of each side and expanding ojr, gives 

d 


lnpi(i/i) = \nuR (Et) - 

-lnZ' + ..., 


dE. 


Incu/j (Et) 


(25) 


where Z' is a constant. The higher-order terms that are indicated by the dots in Eq. (25) 


are proportional to the ratio of Ei to the much larger values of Er, and in the limit of an 
inhnite reservoir, they vanish. Although this calculation is based on finding the mode of 
the probability distribution, rather than the mean, the assumption of an inhnite reservoir is 
sufficient to make the mean and mode agree. Interactions between hnite subsystems will be 
discussed elsewhere, but the effects are proportional to the inverse of the number of particles, 
and are not measurable for macroscopic subsystems. 

Dehning 


/9 = 


d 


dE. 


In Ur (Et) , 


Eq. (25) becomes 


lnpi(z/i) = IncjR {Et) - ( 3Ei{ui) -InZ', 


or 


(26) 

(27) 

(28) 


Pi)"]) = ^WH(BT)exp(-/?Bi(i^,)) . 

Since ur{Et) does not depend on i/i, we can simplify the expression by dehning a new 
constant Z = Z'Iu 2 (^ 1 ^ 2 ), which gives 


PiWi) = -^exp(-/3Ei(vi)) 


(29) 


This is, of course, just the canonical distribution for subsystem 1 for an inverse temperature 


/3 = l/ksT. The partition function Z = Z{l3) in Eq. (29) is given by the normalization 
condition, 

Z = ^exp(-/?Ei(z/i)). (30) 


El 


Note that the only property of the reservoir that enters into this equation is the inverse 
temperature f3, which is the only property of the reservoir needed to determine the thermo¬ 
dynamics of subsystem 1. 
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The average energy of subsystem 1, which can be identified with the thermodynamic 
energy, f/i, is 


um = 5 : Eii.i'i) pii.i'i) 




Z{f3) 


Ei{ui) exp {-(]Ei{ui)). 


(31) 


III.2. Equilibrium between finite systems 


An important concern for practical applications is how the values of for subsystem 


1 in Eq.(29) are affected by replacing the infinite reservoir R, used in Eq. (24), by a hnite 


subsystem 2. The higher-order terms in Eq. (25) no longer vanish, but they are small 


It is not really necessary for Eq.(29) to be valid for all values of Ei{i'\). Because the 
values of Pi(z^i) will only be signihcantly different from zero within the thermal fluctuations 


6 El, it is only necessary for Eq.(29) to be valid over a range of energies of the order of 5Ei. 


As long as the ratio 6 E 1 /U 2 is small, Eq.(29) will be accurate. From a consideration of the 


magnitude of the fluctuations, the condition of validity is that 

6Ei 


Ni 


« 1 . 


(32) 


U2 V «' 2 (A'i + N2) 

This condition will be easily satisfied for N 2 » Ni. For the opposite case of subsystem 1 


being in equilibrium with a relatively small subsystem 2, so that N 2 « N 2 , Eq. (32) becomes 


6 E 1 /U 2 ~ 1 /\/N 2 . This condition can also be satished if subsystem 2 is macroscopic. 

The conclusion is that it doesn’t matter much for the distribution of energies whether 
a subsystem is in equilibrium with a large subsystem or a small one - as long as they 
are macroscopic. This is consistent with general experience of measured temperatures and 
experimental conhrmations of the zeroth law of thermodynamics. 


The next step is to explore the consequences of Eqs. (29), (30), and (31) for the entropy 
through the use of a Massieu function. 


III.3. Massieu functions 

It is well known that the canonical partition function is related to the Helmholtz free 
energy, E = U — TS, hj the equation 

\nZ{P,V,N) = -pE{T,V,N), (33) 
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where T is the temperature, and f3 = l/ksT. The Helmholtz free energy is, of course, the 
Legendre transform of the fundamental relation U = U{S, V, N) with respect to temperature, 
which we will denote as F{T, V, N) = U[T], indicating the Legendre transform by the square 
brackets around the new variable T |26l EZj • 

Although the use of the thermodynamic potentials derived from the fundamental equation 
U = U{S, V, N) through Legendre transforms are familiar to all students of thermodynamics, 
they are not appropriate for calculating the properties of non-monotonic entropy functions. 
The reason is that if S' = S{U, V, N) is not monotonic in U, the function cannot be inverted 
to hnd U = U{S,V, N). However, we can still use Massieu functions, which are Legendre 
transforms of S = S{U, V, iV)[2ni EZ]- It will be particularly useful to dehne a dimensionless 
entropy, S = S/kB, in forming Massieu functions. 

From the differential form of the fundamental relation for dS, we can see that 


(34) 


dS = pdU + pPdV - 0)1 dN, 


where P is the pressure, V is the volume, /i is the chemical potential, and N is the number 
of particles. The inverse temperature f3 is found from the usual equation, which can be 
written as 



(35) 


The Legendre transform (Massieu function) is given by 


^[/3] = S - (3U = -(3{U -TS) = -(3F, 


(36) 


SO that 


S|/?]=liiZ(/?,ViV). 


(37) 


The differential form of the Massieu function S')/?] is then 


dS[(3] = -Udl3 + IdPdV - (df^dN. 


(38) 


This immediately gives us 



(39) 


where the last equality is a well-known thermodynamic identity [26] Wi\. 
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III.4. Inverse Legendre transform of 5[/3] to find S{U) 


To carry out the inverse Legendre transform of S'[/9] to find S{U), use Eq. (39) to find 
U = U (/3). Since ?7 is a monotonic function of /3, even for a non-monotonic density of states, 


we can invert this equation to obtain (3 = (3{U). From Eq. (35), we can find 


S = S[(3] + (3{U)U. 


(40) 


Finally, the entropy with the usual dimensions is given by 

^ = ksS. 


(41) 


III.5. Alternative calculation of the entropy S{U) from d = P{U) 


Since (3 = (3{U) was found by inverting U = U{(3), the entropy can also be calculated by 


integrating the thermodynamic identity in Eq. (35) to find 

'■U 


S = k 


B 


l3{U')dU' 


(42) 


Urnin 

where Umin is the minimum value of the thermodynamic energy U. Both methods produce 
the same results. 

The next section discusses the application of these methods to calculate S{U) for either 
a monotonic entropy or a non-monotonic entropy. 


IV. EXAMPLES OF THE ENTROPY OF QUANTUM SYSTEMS 

This section contains explicit calculations of the entropy of quantum systems, as illus¬ 
trations of the methods described in the previous section. The first example is a system of 
quantum simple harmonic oscillators, which has a monotonic dependence of the entropy as 
a function of energy. . The second example is a system composed of two-level objects. This 
system has a non-monotonic entropy and illustrates how negative temperatures can arise in 
a quantum system. For completeness, I will discuss two kinds of two-level systems. 

A clear distinction should be made between macroscopic systems composed of microscopic 
objects (simple harmonic oscillators, two-level objects, etc.), as discussed in this section, and 
the composite systems discussed in Sections [IT| and [ITU which were composed of macroscopic 
subsystems. 
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IV. 1. A system composed of quantum simple harmonic oscillators 


Consider a system composed of N simple harmonic oscillators. For simplicity, we will 
assume that the frequencies of all oscillators are the same, so that the Hamiltonian is 

HsHO,k = hw{nk + 1/2), (43) 

where h is Planck’s constant, u is the angular frequency, and = 0,1, 2,... The partition 
function of the complete system of N oscillators is well known to be 

ZsHo = exp{-/3hwN/2) (1 - exp{-/3hw))~^ , (44) 


so that 


‘S'sHO [/ 3 ] — In .^SHO 


= --N/3hw - Nln{l - exp{-/3hw)) 


(45) 


The energy is given by the negative partial derivative of ^sHo[/d] with respect to ( 3 , as in 

UsHo = -jhxjN + NhxJ (exp(/3ha;) — 1)”^ (46) 


Eq. d^. 


IV. 1.1. Finding Ssho by Legendre transform 


To make the following equations more compact, dehne 

^SHO = UsHO — -hwV. 


(47) 


The next step is to invert Hsho = ^sho(/^)) to hnd (3 = / 3 {Usro)- To make the notation 
still more compact, we will express the results in terms of a dimensionless energy variable 

f^SHO — 


X = Usno/Nhw = 


Nhu 


which gives 


^ In [1/x + 1] = ^ In 
nw nuj 


1 + a; 


X 


(48) 


(49) 


The inverse Legendre transform, ^sho[/ 9] as a function of x is 


*S'sHo[/^] = --iVln [l/x + 1] + iVln (1 + x)) 


(50) 
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Using Eq. (40), the dimensionless entropy is found to be 


‘S'sHO = fcfi'S'sHo = Nks [-x In a; + (1 + x) ln(l + x)] (51) 

It is readily confirmed that the same result is obtained by numerical integration of the inverse 
temperature. 

In the limit of /3 —>■ oo (or T —)■ 0) the energy goes to its minimum value, f/sno —t huN/2, 
and X —)■ 0 . In this limit, S'sho —t 0 , as expected from the third law of thermodynamics 
(Nernst theorem). 


IV.2. A system of two-level quantum objects 


The next example illustrates the thermodynamics of a system with a bounded energy 
spectrum and a non-monotonic entropy. This system displays negative temperatures. 
Consider a collection of N two-level quantum objects, 

N 

.^2-level ^ ^ 
k=l 

where e is the energy difference between the two levels in each object, and = 0 or 1. The 
partition function Z 2 -ievei can be found by standard methods. 

^ 2 -ievei = [1 + exp(-/3e)]^ (53) 


This gives 


and, using Eq. (39), 


S2-ievei|/3| = 111 Zj-fevei = A' Ih [l + exp(-^e)l, 
= A'e (exp(/?6) + 1)"‘. 


(54) 

(55) 


IV.2.1. Finding S 2 -ievei{U) by Legendre transform 


As before, we simplify the notation by dehning a dimensionless energy 

fA-level 


y = 


Ne 


(56) 


The same type of calculation used for the simple harmonic oscillators gives the inverse 
temperature as 

1 n 1 1 n - - 7/1 

(57) 


1 , 

■f 

1 , 

1 - 

1 

1 — 1 

^n 

- - 1 

= - In 

e 

[y J 

e 

L y \ 
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A little algebra then gives 


*S' 2 -ievei = -Nks [y\ny+ {1- y) ln(l - y)] (58) 

As was the case for the system composed of simple harmonic oscillators, integration of the 
inverse temperature produces the same result as the Legendre transform. 

*S' 2 -ievei has positive temperatures for y < 0.5, and negative temperatures for y > 0.5, as 
expected. The entropy is symmetric for y 1 — y. The entropy goes to zero in the limits 
?/ —;• 0 and 1 / —t 1, also as expected. 


IV.3. Independent Ising spins 


Since a system composed of non-interacting Ising spins is isomorphic to the two-level 
system discussed above, its entropy can be obtained with little effort. 

Consider the Hamiltonian 

N 

hfspins ^ ^ ^ j 1 

where Uj takes on the values -|-1 and —1, and the parameter b represents an applied magnetic 
held. The mapping between the variable Uj, which take on the values 0 and 1, and the spins 

is 

CTj = 2nj — 1, (60) 


and the energies map as 


h^2-level Cgpins 1 


Ne 2Nb ^ 2' 


Dehning a dimensionless energy for the spin system by 




Eq. (61) becomes 


r = 


y = 


spins 


Nb ’ 

r -|- 1 


(61) 


(62) 


(63) 


where y is the dimensionless energy for the two-level system dehned in Eq. (56). This gives 
the entropy of the spin system as 

Nks 


*S*spins 


[(1 — r) ln(l — r) 


-I- (1 -I- r) ln(l -I- r) — 2 In 2] . 


(64) 
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S'spins has positive temperatures for r < 0, and negative temperatures for r > 0, as 
expected, with an obvious symmetry for r ■H- —r. The entropy goes to zero in the limits 
r —)■ — 1 and r —)■ +1, also as expected. 

A slight modification of the derivations given above would also allow the calculation of 
the entropy of a system composed of two-level systems or simple harmonic oscillators with 
arbitrary distributions of energy-level spacings. 

Now that these examples of derivations of the entropy of quantum systems are complete, 
they can be tested against exact results and compared with the predictions of the volume 
entropy advocated by opponents of the concept of negative temperatures. 


V. COMPARISON OF 52-ievei WITH THE VOLUME ENTROPY 


Comparisons with the volume entropy as used by its advocates are somewhat ambiguous 
because they present their ideas in terms of a non-analytic expression for the entropy of a 
quantum system[Tl H]. However, their basic claim is that the logarithm of the integral of 
the density of states gives the true entropy, which always has positive temperatures. The 
volume entropy is calculated in the next section using the same procedure used in Ref. HI 
The predictions of the volume entropy are then tested against exact results and compared 


with the expressions of the entropy derived in Section IV 


V.l. Calculation of the volume entropy for a system composed of two-level objects 


In the notation of Ref. HI 52_ievei is related to the density of states u{y) by 

*S2-ievei = lni:n(|/) 

= -N[y\ny+{l-y) ln(l - y)] , 


(65) 


where the explicit dependence on the dimensionless energy y is taken from Eq. (58). Since 


the dimensionless entropy is being expressed in terms of a dimensionless energy variable, 
y = U2-\eve\l Ne, there is no need for the energy parameter that Ref. HI included in the 
logarithm. 

The dimensionless volume entropy is given by 


Sciy) = lnH(2/), 


( 66 ) 
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where 


n(!/) = / ^iy' 




(67) 


Carrying out the integral in Eq. (67) (numerically) gives Sq- After this is done, the fact 
that 

dn{y) 


^iy) = 


dy 


( 68 ) 


gives a simple expression for the inverse temperature. 

These explicit expressions for the volume entropy and its associated temperature are 
tested against exact results in the next subsection. 

V.2. Comparison between predictions of the partition of energy between subsys¬ 
tems in equilibrium for different entropies 

Consider thermal equilibrium between two two-level subsystems with the same level spac¬ 
ing e, but differing in size by a factor of /. Subsystem 1 contains fN two-level objects, and 
subsystem 2 is 1// smaller, with N objects. Even without calculations, it is obvious that 
equipartition of energy between the two subsystems in equilibrium requires 


fN 


U2 

N' 


(70) 


The prediction of S' 2 _ievei comes from setting the temperatures equal. Using Eq. (57), 

fjNe 


Pj = - In 


_ f72-level j' 


- 1 


1 

r 1 1 

= - In 

— 1 

e 

m \ 


(71) 


for j = 1 or j = 2, and /i = /, while /2 = 1. Since e is the same in both subsystems, it is clear 


that the condition of equilibrium is yi = ?/ 2 , which is equivalent to Eq. (70). This conhrms 


that S' 2 -ievei correctly predicts equipartition for these subsystems with non-monotonic oj{y). 

The numerical computations for the equilibrium predictions of the volume entropy have 
been have been carried out for test cases and plotted in Fig. For comparison, the exact 


condition for equilibrium, which is a straight line given by Eq. (70), is also shown. This 


exact equilibrium condition is identical to that predicted by S' 2 -ievei, as derived in Section 
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IV.2 For values of the energy less than half the maximum energy, all predictions are in 


good agreement with the equipartition of energy. However, for values of the energy more 
than half the maximum energy, the predictions of the volume entropy deviate signihcantly 
from the correct equilibrium conditions. 

The origin of the errors made by the volume entropy can be traced back to the size 
dependence of the temperature Tq given by the volume entropy: the larger the subsystem, 
the higher Tg- Since subsystem 1 is larger than subsystem 2, it also has a higher value of 
Tg at equilibrium than subsystem 2. Equal values of Tg in the two subsystems can only be 
achieved if subsystem 2 has a higher energy than appropriate for equilibrium, leading to the 
erroneous predictions of the volume entropy shown in Fig. 



FIG. 1. Plot of the predictions of the volume entropy are compared with those of the entropy given 


by 5 ' 2 -level in Eq. (58) for the distribution of energy between two subsystems of two-level objects 
in thermal equilibrium. The exact condition of equipartition of energy is indicated by the straight 
line, which is also the prediction of 52_ievei- The three curves giving the volume entropy predictions 
are for subsystem sizes [A^i, A^ 2 ] of [10,1], [100,10], and [1000,100]. In each case, the value of N 2 is 
given in the legend. Larger deviations of the curves for the volume entropy from the exact results 
correspond to larger subsystems. Increasing the ratio of the subsystem sizes also increases the 
violation of equipartition of energy. 
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VI. CONCLUSIONS 


I have demonstrated that the entropy of a macroscopic quantum system is a continuous 
function of the thermodynamic energy, as opposed to the step functions and delta functions 
proposed by other authors [H HJ |6l IMl 1 ^ - 


The expressions for the entropy of simple harmonic oscillators in Eq. (51) and two-level 


systems in Eqs. (58) and (64) satisfy the postulates of thermodynamics. They are funda¬ 
mental thermodynamic relations that completely characterize the correct thermodynamic 
properties of these systems. Because the fundamental relation of a thermodynamic system 
can be shown to be unique, any valid dehnition of the entropy of macroscopic quantum 
systems must be equivalent to the one presented in this paper. 

An interesting feature of the dehnition of entropy presented here is that subsystems with 
incommensurate energy eigenvalues can be in equilibrium with each other. It is, of course, 
well known experimentally that differences in microscopic energy-level spacings do not pre¬ 
vent equilibration, but this fact has not been obvious from earlier proposed expressions for 
the entropy of macroscopic systems. 

Finally, I have shown that the entropy of two-level objects conhrms the validity of the 
concept of negative temperatures. 
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